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ABSTRACT 

The magnetorotational instability is thought to be responsible for the generation of 
magnetohydrodynamic turbulence that leads to enhanced outward angular momentum 
transport in accretion discs. Here, we present the first formal analytical proof showing 
that, during the exponential growth of the instability, the mean (averaged over the 
disc scale-height) Reynolds stress is always positive, the mean Maxwell stress is always 
negative, and hence the mean total stress is positive and leads to a net outward flux of 
angular momentum. More importantly, we show that the ratio of the Maxwell to the 
Reynolds stresses during the late times of the exponential growth of the instability 
is determined only by the local shear and does not depend on the initial spectrum 
of perturbations or the strength of the seed magnetic. Even though we derived these 
properties of the stress tensors for the exponential growth of the instability in incom- 
pressible flows, numerical simulations of shearing boxes show that this characteristic 
is qualitatively preserved under more general conditions, even during the saturated 
turbulent state generated by the instability. 
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1 INTRODUCTION 

Magnetohydrodynamic (MHD) turbulence has long 
been considered responsible for angular momentum trans- 
port in accretion discs su rrounding astrophysical objects 
iShakura fc Sunvaevl Il973l) . Strong support for the im- 
portance of magne tic fields in accre t ion d iscs followed 
the realization by iBalbus fc Hawle'vl (Il99ll) that lami- 
nar flows with radially decreasing angular velocity pro- 
files, that are hydrodynamically stable, turn unstable 
when threaded by a weak magnetic field. Since the dis- 
covery of this magne torotational instability (MRI), a 



is best appreciated by examining the equation for the dy- 
namical evolution of the mean s pecific angular momen- 
tum of a fluid element (se e, e.g.. IBalbus fc Hawlevl Il998l : 
IBalbus fc Papaloizoul Il999h . Defining this quantity for a 
fluid element with mean density p as 1 = rpv^ we can write 
in cylindrical coordinates 



d t l+V-(lv) 



— V-r{B 4 ,B)+ — V-r5B 4> 5B-V-r 

47T 47T 



(1) 



varie ty of local (e.g., |HawJgj^_GarmTrte J- J^^albjisJ W9 



1 99(|_jBnindenbure^fOT^ 

iBrandenburd bOOljiJSano^jnutsuka. Tnme£ J _fc_Sj^on6 
and global (e.g.. lHawlevll2000l. l200ll: IStone fc Pringl 



1995; 



2004) 



2001) 



numerical simulations have shown that its non-linear evolu- 
tion gives rise to a turbulent state characterized by enhanced 
Reynolds and Maxwell stresses, which in turn lead to out- 
ward angular momentum transport. 

The relevance of the Reynolds and Maxwell stresses in 
determining the dynamics of a magnetized accretion disc 



where v,B, 8v, and SB stand for the means and fluctua- 
tions of the velocity and magnetic fields, respectively, and 
the bars denote suitable averages. It is clear that the pres- 
ence of mean magnetic fields or of non-vanishing correla- 
tions between the fluctuations in the magnetic or velocity 
field can potentially allow for the specific angular momen- 
tum of a fluid element to change. In the absence of strong 
large-scale magnetic fields, the last two terms on the right 
hand side will dominate and we can simplify equation |0 as 



dtl + V-(lv) 



(2) 



where the vector T characterizes the flux of angular mo- 
mentum. Its components are related to the Reynolds and 
Maxwell stresses, Rij = pbvi Svj and Mij = SBi SBj/4-K via 
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It is straightforward then to see that in order for matter in 
the disc to accrete, i.e., to loose angular momentum, the sign 
of the mean total stress, T r< f,, must be positive. Note that in 
the MRI-literature the Maxwell stress is defined as the neg- 
ative of the correlations between magnetic field fluctuations. 
We have chosen instead to use a definition in which the sum 
of the diagonal terms has the same sign as the magnetic 
energy density. 



In a differentially rotating MHD turbulent flow, 
the sign of the stress component T r( f, will depend on 
the mechanism driving the turbulence (presumably 
the MRI), the mechanism mediating the energy cas- 
cade between different scales, and the dissipative pro- 
cesses that lead to satur ation. Nevertheles s , mec h anica l 
analogies of the MRI jBalbus fc Hawlevl Il99l 11998 



KatorFukue' fc Mineshigelll99Sl: iBrandenburg fc Campbelll 
19971 : lBalbusll2003i) as well as anal yses involving the excita- 
tion of single- wavenumber modes llBaIJ>us_^JIawleyj |l992|, 
2002 ; iNaravan. Quataert. Igumenshchev. fc Abramowiczl 
2002, see also §4) suggest that angular momentum in 
transported outwards even during the linear phase of the 
instability. 



In this paper, we derive analytic expressions that re- 
late the dynamics of the MRI-driven fluctuations in Fourier 
space with the mean values, averaged over the disc scale- 
height, of the different stresses in physical space. This allows 
us to provide the first formal analytical proof showing that 
the MRI leads to a positive mean total stress T r( j, which, 
in turn, leads to a net outward transport of angular mo- 
mentum. Within this formalism, we demonstrate that the 
Reynolds stress is always positive, the Maxwell stress is al- 
ways negative, and that the absolute value of the Maxwell 
stress is larger than the Reynolds stress as long as the flow 
in the absence of magnetic fields is Rayleigh-stable. More- 
over, we uncover a robust relationship between the Maxwell 
and Reynolds stresses as well as between the magnetic and 
kinetic energy of the MHD fluctuations during the late times 
of the linear phase of the instability. Specifically, we show 
that the ratio of the (absolute value of the) Maxwell and 
Reynolds stresses is equal to the ratio of magnetic and ki- 
netic energy densities, and that both ratios depend only on 
the value of the local shear characterizing the flow. 



The rest of the paper is organized as follows. In §2 we 
state our assumptions. In §3 we present the complete solu- 
tion to the MRI-eigenvector problem in Fourier space. We 
pay particular attention to the characterization of the com- 
plex nature of the MRI-eigenvectors. In §4, we present the 
formalism to derive analytic expressions for the mean phys- 
ical Reynolds and Maxwell stresses in terms of the fluctu- 
ations in spectral space. We show that the MRI leads to a 
net outward transport of angular momentum. We also study 
there the different properties of the MRI-driven stresses 
during the late times of their exponential growth. In §5 
we compare these stress properties with similar properties 
found in previous numerical simulations that addressed the 
non-linear turbulent regime using shearing boxes. We also 
present there our conclusions. 



2 ASSUMPTIONS 

This paper is concerned with the signature of the ax- 
isymmetric MRI in the mean values (averaged over the 
disc scale-height) of the Reynolds and Maxwell stress ten- 
sors. In particular, we consider a cylindrical, incompress- 
ible background, characterized by an angular velocity pro- 
file fl — Q(r)z, threaded by a weak vertical magnetic field 
B — B z z. In order to address this issue, we work in the 
shearing sheet approximation, which has proven useful to 
understand the physics of disc phenomena when the scales 
involved are smaller than the disc radius. 

The shearing sheet approximation consist of a first order 
expansion in the variable r — ro of all the quantities charac- 
terizing the flow at the fiducial radius ro. The goal of this ex- 
pansion is to retain the most important terms governing the 
dynamics of the MHD fluid in a locally-Cartesian coordinate 
system co-orbiting and corrotating with the background flow 
with local (Eulerian) velocity v = ro f2o</>- (For a more de- 
tailed di scussion regarding the shearing sheet approxima- 
tion, see lGoodman fc Xul fl994 and references therein.) 

The equations for an incompressible MHD flow in the 
shearing sheet limit are given by 



dv 

Tit 



+ {v-V)v = -2fl xv + gfioV(r - r ) 2 

l„f n B 2 \ (B-V)B 
--V PH + — 

p \ 8ty J 4np 



dB 
~dt 



+ (v-V)B = (B-V)v 



(4) 
(5) 



where P is the pressure, p is the (constant) density, the 
factor q = — dlnf2/dlnrj ro parametrizes the magnitude of 
the local shear, and we have defined the (locally-Cartesian) 
differential operator 



,9 d> d _ d 
or ro U(f> Oz 



(6) 



where f, </>, and z are, coordinate- independent, orthonormal 
vectors corrotating with the background flow at ro. 

In what follows, we focus our attention on fluctuations 
that depend only on the vertical coordinate. These types 
of fluctuations are known to h ave the fastest growth rates 
JBalbus fc Hawlevl If992l . Il998l) and will, therefore, consti- 
tute the most important contributions to the correspond- 
ing Reynolds and Maxwell stresses during the exponential 
growth of the instability. The equations governing the dy- 
namics of these fluctuations can be obtained by noting that 
the velocity and magnetic fields given by 

v = Sv r (z)r + [-q£l (r — ro) + 5v<p(z)]4> + Sv z (z)z , (7) 
B = 5B r (z)r + 8B^(z)^ + [B z +SB z (z)]z, (8) 

constitute a family of exact, non-linear, solutions to the one- 
dimensional incompressi ble MHD equation s in th e shearing 
sheet limit. As noted in iGoodman fc Xul Jl994), the only 
non-linear terms, which are present through the perturbed 
magnetic energy density, are irrelevant in the incompressible 
case under study (i.e., the total pressure can be found a 
posteriori using the condition V-u = 0). 

Due to the divergenceless nature of the disturbances 
under consideration, the fluctuations in the vertical coordi- 
nate, 5vz(z) and 5B z (z), reduce to a constant. Without loss 
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of generality, we take both constants to be zero . We can 
further simplify the system of equations @ an d @ by re- 
moving the background shear flow i! s hcar = — q£lo(r — ro)(j>. 
We obtain the following set of equations for the fluctuations 



9 x 
dt 5Vr 

m Sv * 

d_ 

dt 
o_ 

dt 



8B r = 



B 
2Q.o5va, + 7— SB r , 
Aitp oz 

-(2-q)Q, 5v r +^^-8B 4> , 
<mp oz 

Bz-7T-8v r , 

oz 

- d 

—qQoSBr + B z — 8vd, , 
oz 



(9) 
(10) 

(11) 
(12) 



where the first term on the right hand side of equation Hot 
is related to the epicyclic frequency k = \/2(2 — q) f2o , at 
which the flow variables oscillate in an unmagnetized disc. 

It is convenient to define the new variables Sbi = 
SBi I ytiirp for i = r, <f>, and introduce dimensionless quanti- 
ties by considering the characteristic time- and length-scales 
set by 1/fio and B z /(^/47rp flo). The equations satisfied by 
the dimensionless fluctuations, 8vi, Sbi, are then given by 



df5l!r 

df8b r 
djSbj, 



2SV4, + diSb r , 

— (2 — q)8v r + dgSbj, . 

dz5v r , 

—q5b r + diSv^ , 



(13) 
(14) 
(15) 
(16) 



where t and z denote the dimensionless time and vertical 
coordinate, respectively. 

In order to simplify the notation, we drop hereafter the 
tilde denoting the dimensionless quantities. In the rest of the 
paper, all the variables are to be regarded as dimensionless, 
unless otherwise specified. 



3 THE EIGENVALUE PROBLEM FOR THE 
MRI : A FORMAL SOLUTION 

In this section we provide a complete solution to the set 
of equations 113H - 116t in Fourier space. Taking the Fourier 
transform of this set with respect to the z-coordinate, we 
obtain the matrix equation 



d t S(k n ,t) = L5(k„,t) , 

where the vector S(k n ,t) stands for 



5(kn,t) = 



Sv r (k n , t) 
8v${k n ,t) 
Sb r (k n ,t) 

Sb$(k n ,t) . 
and L represents the matrix 



(17) 



(18) 






2 


ik n 





(2-9) 








ik n 


















-Q 






(19) 



The functions denoted by f(k n , t) correspond to the Fourier 
transform of the real functions, f(z, t), and are defined via 



1 Note that the fluctuations studied in iGoodman fc Xul dl9& 
are a particular case of the more general solutions J7J and 181 , 



f(kn,t) 



1 

2H 



f(z,t)e lknZ dz . 



(20) 



where H is the (dimensionless) scale-height and k n is the 
wavenumber in the z-coordinate, 



_ rnv 
fc„ = — . 



(21) 



with n being an integer number. Here, we have assumed 
periodic boundary conditions at z — ±H . 

In order to solve the matrix equation \Y7} . it is conve- 
nient to find the base of eigenvectors, {ej } with j = 1, 2, 3, 4, 
in which L is diagonal. This basis exists for all values of the 
wavenumber k n (i.e., the rank of the matrix L is equal to 4, 
the dimension of the complex space) except for k n = and 
k n = %/2~9- l n thi s base, the action of L over the set {e-, } is 
equivalent to a scalar multiplication, i.e., 

idiag Bj = ffj ej for j = 1,2,3,4, (22) 

where {oj} are complex scalars. 

3.1 Eigenvectors 

In the base of eigenvectors, the matrix L has a diagonal 
representation Ldiag = diag(cri, 02, 03, 04). The eigenvalues 
{<Tj}, with j = 1,2,3,4, are the roots of the characteristic 
polynomial associated wi th L, i.e., the dispersion rela tion 
associated with the MRI fealbus fc Hawlevlll99illl998l) . 



(kl + a]) 2 + 2(2 ~q){kl + a]) 
and are given by 

1/2 



4A£ = 0, 



(23) 



= ± -A± VA 



(24) 



where we have defined the quantities A and A such that 



kl + (2 - q) , 



(25) 

A ee (2-g) 2 + 4fc^. (26) 

For the modes with wavenumbers smaller than 

k BH = ^, (27) 

i.e., the largest unstable wavenumber for the the MRI, the 
difference vA — A is positive and we can define the "growth 
rate" 7 and the "oscillation frequency" uj by 



(Va-a 

v / A + A 



1/2 
1/2 



(28) 
(29) 



both of which are real and positive (for all positive values 
of the parameter q). This shows that two of the solutions of 
equation (12M> are real and the other two are imaginary. We 
can thus write the four eigenvalues in compact notation as 



<T2 



03 



and 



0-4 



(30) 



The set of normalized eigenvectors, {e CTj }, associated 
with these eigenvalues can be written as 



for 3 = 1,2,3,4, 



(31) 



where 
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/(&n) — 



{kl + <rf)/2 
-i2k„<Tj/(kn + a]) j 



(32) 



and the norms are given by 
r 4 1 1/2 

ii e j|i = E< 



(33) 



where e^ is the Z-th component of the (unnormalized) eigen- 
vector associated with the eigenvalue o~j. This set of four 
eigenvectors {e CT .}, together with the set of scalars {<Jj}, 
constitute the full solution to the eigenvalue problem de- 
fined by the MRI. 

The roots Oj of the characteristic polynomial are not 
degenerate. Because of this, the set {e CT . } constitutes a basis 
set of four independent (complex) vectors that are able to 
span C 4 , i.e., the space of tetra-dimensional complex vectors, 
for each value of k n , provided that 7 and u) are given by 
equations I28H and 1291 . respectively. Note, however, that 
they will not in general be orthogonal, i.e., e CTj . ■ e CT ., 7^ for 
j 7^ j' . If desired, an orthogonal basis of eigenvectors can 
be constructed usin g the Gram-Schmidt ort hogonalization 
procedure (see, e.g.. iHoffman fc Kunzdll97ll) . 



3.2 Properties of the Eigenvectors 

Despite the complicated functional dependence of the 
various complex eigenvector components on the wavenum- 
ber, k n , some simple and useful relations hold for the most 
relevant (unstable) eigenvector. Figure^shows the four com- 
ponents of the eigenvector e 7 as a function of the wavenum- 
ber for a Keplerian profile (q = 3/2) and illustrates the fact 
that the modulus of the different components satisfy two 
simple inequalities 



> 
> 



(34) 
(35) 



for all values of < k n < feBH- These inequalities do indeed 
hold for all values of the shear parameter < q < 2. 

We also note that equation 1321 1 exposes a relationship 
among the components of any given eigenvector e CT . . It is 
immediate to see that 

— el. Oik^ 

for j= 1,2,3,4. (36) 



ei j el j kl + O j 

In particular, the following equalities hold for the compo- 
nents of the unstable eigenvector, e 7 , at the wavenumber 



fcmax = — W - — 1 , (37) 

for which the growth rate is maximum, 7 max = ?/2, 

(38) 



3-j,(/Cmax) — G-)/(&max) — ^ 



(ha 



(39) 



As we show in the next section, the inequalities l|34|l 
and 13511 . together with equations 1381 and 13911 . play an 
important role in establishing the relative magnitude of the 
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Figure 1. The components of the normalized unstable eigenvec- 
tor e 7 defined in equation 1311 . The vertical dotted lines denote 
the wavenumber corresponding to the most unstable mode, fc max , 
(eq. 1371 ), and the largest unstable wavenumber, fceHi ( e Q- H3 )■ 



different mean stress components and mean magnetic en- 
ergies associated with the fluctuations in the velocity and 
magnetic fields (see also Appendix lB"t. 

Finally, we stress here that the phase differences among 
the different eigenvector components cannot be eliminated 
by a linear (real or complex) transformation. In other words, 
it is not possible to obtain a set of four real (or purely imag- 
inary, for that matter) linearly independent set of eigenvec- 
tors that will also form a basis in which to expand the gen- 
eral solution to equation 1171 . Taking into consideration the 
complex nature of the eigenvalue problem in MRI is crucial 
when writing the physical solutions for the spatio-temporal 
evolution of the velocity and magnetic field fluctuations in 
terms of complex eigenvectors. As we discuss in the next 
section, this in turn has a direct implication for the expres- 
sions that are needed in the calculation of the mean stresses 
in physical (as opposed to spectral) space. 

3.3 Temporal Evolution 

We have now all the elements to solve equation 1171 . In 
the base defined by {e CT;j }, any given vector S(k n ,t) can be 
written as 

4 

S(k n ,t) = ^2a j [k n ,t)e ai , (40) 
3=1 

where the coefficients a,j(k n ,t), i.e., the coordinates of 
S(k n ,t) in the eigenvector basis, are the components of the 
vector a(k n ,t) obtained from the transformation 

a(k n ,t) = Q- 1 6(k n ,t) . (41) 

The matrix Q~ x is the matrix for the change of coordinates 
from the standard basis to the normalized eigenvector basis 
and can be obtained by calculating the inverse of the matrix 

Q = [e CT1 e CT2 e CT3 e CT J . (42) 

Multiplying equation (11711 at the left side by Q _1 and 
using the fact that 

Ldiag^Q-'LQ, (43) 
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we obtain a matrix equation for the vector a(k n ,t), 
dta(k„,t) = Ldiag a(k n ,t) , 


(44) 


mergies 


0.5 
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which can be written in components as 




<« 


0.4 






-M r0 


dtOij(k n ,i) = aj aj(k n ,t) with j = 1,2,3,4. 
The solution of these equations is then given by 


(45) 


Stresses 


0.3 
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R r* 
•9 

E u 


aAk* t) — aAk„ 0)e' Tjt with 1 = 12 3 4 
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CD 
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We can finally write the solution to equation (11 71 


as 


zed p 




i / / 
1 / / 


\ \ 




4 

S(k n ,t) = 2J aj(fc„,0) e CTj ' e CTj . , 

3=1 


(47) 


Normal! 


0.1 


-i / / 

' ' / 
- / / 

s^, 1 







where {<Jj} and {e CT .}, for j — 1,2,3,4, are given by equa- 
tions (13011 and 1311 . and the initial conditions a(k n ,0) are 
related to the initial spectrum of fluctuations, S(k n ,0), via 
a(k n ,0)=Q- 1 5{k n ,0). 



4 NET ANGULAR MOMENTUM 
TRANSPORT BY THE MRI 

Having obtained the solution for the temporal evolution 
of the velocity and magnetic field fluctuations in Fourier 
space we can now explore the effect of the MRI on the mean 
values of the Reynolds and Maxwell stresses. 



4.1 Definitions 

The average over the disc scale-height, 2H , of the prod- 
uct of any two physical quantities, f(z, t) g(z, t), 



(48) 



fg(t) = I f(z,t)g(t,z)dz : 



can be written in terms of their corresponding Fourier trans- 
forms, f(k n ,t) and g(k n ,t), as (see Appendix IXt 



fg(t) = 2^ Re[ /(*».*) 5* (*»,*)] 



(49) 



Here, Re[ ] stands for the real part of the quantity between 
brackets, the asterisk in g*(k n ,t) denotes the complex con- 
jugate, and we have considered that the functions f(z,t), 
and g(z, i) (both with zero mean) are real and, therefore, 
their Fourier transforms satisfy /(— k n ,i) = f*(k n ,t). 

Using equation 14911 . we can write the mean values of 
the quantities 



Mij(z,t) 



8vi(z,t) Svj(z,t) , 
Sbi(z,t) 8bj(z,t) , 



with i,j = r, (f>, as 

Rv(t) 



2^ Re[Svi{k n ,t)S'v;(k n ,t)'_ 

n=l 

oo 

M l3 {t) = 2J2 Re[8h(k n ,t)Sb*(k n ,t)} 



(50) 
(51) 



(52) 
(53) 
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Figure 2. The per-k contributions, 1Z r ^(k n ), Mnf>(kn), Exikn), 
and £jy-(fc n ), associated with the corresponding mean physical 
stresses 1541 and 1551 and mean physical magnetic and kinetic 
energy densities, denned in Appendix^] eqs. 1B2I and IB1I , The 
vertical dotted lines denote the wavenumber corresponding to the 
most unstable mode, fc ma x, (eq. 1371 ). and the largest unstable 
wavenumber, fceH; ( e 1- 1271 1. 



4.2 Properties of the MRI-Driven Stresses 

At late times, during the exponential growth of the in- 
stability, the branch of unstable modes will dominate the 
growth of the fluctuations and we can write the most impor- 
tant (secular) contribution to the mean stresses by defining 



JVbh 

R r4 ,(t) = 2^n r4 ,{k n )\a 1 \ 2 e 2lt + ... , 

n=l 

JVbh 

M r4> {t) = 2^ Mr4>{k n )\a 1 \ 2 e 2 '> t + ... , 



(54) 
(55) 



where Abh is the index associated with the largest un- 
stable wavenumber (i.e., the mode labelled with highest 
k n < ^eh) 2 , the dots represent terms that grow at most 
as fast as e 7 *, and we have defined the functions 



n 

and 

Mr<j> 



Re[e* e 2 *] 



ei 



Refe 3 e 4 * 
(M_ lleill 2 



(56) 



(57) 



Note that these functions are not the Fourier transforms 
of the Reynolds and Maxwell stresses, but rather represent 
the contribution of the fluctuations at the scale k n to the 
corresponding mean physical stresses. We will refer to these 
quantities as the per-k contributions to the mean. 

The complex nature of the various components of the 
unstable eigenvector, together with the inequalities 11340 and 
1351 . dictate the relative magnitude of the per-k contribu- 
tions associated with the Maxwell and Reynolds stresses and 
the magnetic and kinetic energy densities (see Appendix ITU 



where the temporal evolution of the fluctuations in Fourier 
space, Sv r (k n ,t), 5v^{k n ,t), Sb r (k„,t), and 5b^(k n ,t), is gov- 
erned by equation 1471 . 



2 Extending the summations to include the non-growing modes 
with k n > &bh would only add a negligible oscillatory contribu- 
tion to the mean stresses. 
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for a detailed discussion of the latter case). Figure [5] shows 
the functions lZ r cj>{k„) and M r <j>{kn) for a Keplerian pro- 
file (q = 3/2). It is evident from this figure that, in this 
case, the per-k contribution of the Maxwell stress is always 
larger than the the per-k contribution corresponding to the 
Reynolds stress, i.e., — M r <t>(k n ) > 1Z r ,p(k n ). This is indeed 
true for all values of the shear parameter < q < 2 (see 
below) . 

The coefficients ei for j = 1,2,3,4 in equations I56H 
and 15711 are the components of the (normalized) unstable 
eigenvector given by equation 1321 with o\ — 7. We can 
then write, the mean values of the Reynolds and Maxwell 
stresses, to leading order in time, as 



JVbh 



M rrt ,(t) 



JVbh 



2, |Ol| 2 7 t 

j ||ei|| 2 ' 

ik'i |aj| 2 2t 
e 



ri k l + 7 2 ||ei| 



(58) 



(59) 



Equations 15811 and 159H show explicitly that the mean 
Reynolds and Maxwell stresses will be, respectively, positive 
and negative, 



Rrfilt) > 



and 



M r <t>{t) < . 



(60) 



This, in turn, implies that the mean total MRI-driven stress 
will be always positive, i.e., 



T r<j> (t) = Rr^t) - M r t(t) > , 



(61) 



driving a net outward flux of angular momentum as dis- 
cussed in 

It is not hard to show now that the magnitude of the 
Maxwell stress, —M r ^(t), will always be larger than the 
magnitude of the Reynolds stress, R r( f,(t), provided that the 
shear parameter is q < 2. In order to see that this is the case, 
it is enough to show that the ratio of the per-k contributions 
to the Reynolds and Maxwell stresses, defined in equations 
1561) and 1571 . satisfy 



-Mr4,(k n ) 



4A:r, 



(fc2+ 7 2 



>1. 



(62) 



for all the wavenumbers k n . Adding and subtracting the 
factor (fc 2 +7 2 ) 2 in the numerator and using the dispersion 
relation 1231 we obtain, 



-Mr#(k„) 

TZ-r4> (^n) 



l + 



2(2 - q) 

kl + 7 2 



(63) 



which is clearly larger than unity for all values of k n provided 
that q < 2. It is then evident that the mean Maxwell stress 
will be larger than the mean Reynolds stress as long as the 
flow is Rayleigh-stable, i.e., 



M r $(t) > Rr4,(t) for < q < 2 . 



(64) 



This inequality provides analytical support to the re- 
sults obtained in numerical simulations, i.e., that the 
Maxwell stress constitutes the major contribution to 
the total stress in magnetized accretion discs (see, e.g., 
lHawlev. Gammie. fc Balbuslll99'5l) . 

We conclude this section by calculating the ratio 
— M r <t>(t) / R r( f,{i) at late times during the exponential growth 
of the instability. For times that are long compared to the 
dynamical time-scale, the unstable mode with maximum 



growth dominates the dynamics of the mean stresses and the 
sums over all wavenumbers in equations 1581 and 159H can be 
approximated by a single term corresponding to fc n — fc max . 
In this case, we can use equations 138H and 1391 to write 



lim ■ 

t»i 



-M r<t ,(t) 
R r<i >(t) 



Re[e 7 e 7 *] 
Refe^f*] 



(65) 



This result shows explicitly that the ratio — M r ^(t) / Rr^(t) 
depends only on the shear parameter and not on the mag- 
nitude or even the sign of the magnetic field. These same 
conclusions can be drawn for the ratio between the mag- 
netic and kinetic energies associated with the fluctuations 
(see Appendix [B] for a detailed discussion). 



5 DISCUSSION 

In this paper we have studied the properties of the mean 
Maxwell and Reynolds stresses in a differentially rotating 
flow during the exponential growth of the magnetorotational 
instability and have identified its signature in their temporal 
evolution. In order to achieve this goal, we obtained the 
complex eigenvectors associated with the magnetorotational 
instability and presented the formalism needed to calculate 
the temporal evolution of the mean Maxwell and Reynolds 
stresses in terms of them. 

We showed that, during the phase of exponential 
growth characterizing the instability, the mean values of 
the Reynolds and Maxwell stresses are always positive and 
negative, respectively, i.e., R r< p(t) > and M rt f,(t) < 0. 
This leads, automatically, to a net outward angular mo- 
mentum flux mediated by a total mean positive stress, 
T r4 ,(t) = R rt j,(t) - M rtf ,(t) > 0. We further demonstrated 
that, for a flow that is Rayleigh-stable (i.e., when q = 
— dlnSl/dlnr < 2), the contributions to the total stress asso- 
ciated with the correlated magnetic fluctuations are always 
larger than the contributions due to the correlations in the 
velocity fluctuations, i.e., —M rr f,(t) > R r<j> {t). 

We also proved that, during the late times of the lin- 
ear phase of the instability, the ratio of the Maxwell to the 
Reynolds stresses simply becomes 

lim -ft*/') = tzl . (66) 

This is a remarkable result, because it does not depend on 
the initial spectrum of fluctuations or the value of the seed 
magnetic field. It is, therefore, plausible that, even in the 
saturated state of the instability, when fully developed tur- 
bulence is present, the ratio of the Maxwell to the Reynolds 
stresses has also a very weak dependence on the properties of 
the turbulence and is determined mainly by the local shear. 

For shearing box simulations with a Keplerian ve- 
locity profile, the ratio of the Maxwell to the Reynolds 
stresses in the saturated state has been often quoted to 
be constant indeed (of order ~ 4), almost independent 
of the setup of the simulation, the initial conditions, 
and the boundary conditions. This is shown in Figure |5] 
where we plot the correlation between the Maxwell stress, 
M r <t>, and the Reynolds stress, R r ^>, at saturation, for 
a number of numerical simulations of shearing boxes, 
with Keplerian velocity profil es but different initial con- 
ditions (data points are from lHawlev. Gammie. fc Balbusl 
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Figure 3. The correlation between the Maxwell stress, M r $, and 
the Reynolds stress, -R r ^, at saturation obtained in numerical sim- 
ulations of shearing boxes with Keplerian velocity profiles but dif- 
ferent initial conditions. Filled circles correspond to simulations 
with zero net vertical magnetic flux. Open circles correspond to 
simulations with finite net vertical magnetic flux. Squares corre- 
spond to simulations with an initial toroidal field. (See the text 
for references.) The dashed line is the correlation often quoted in 
the literature. The solid line corresponds to the late-time ratio 
of the two stresses during the exponential growth of the MRI, as 
predicted by equation 1661 for q = 3/2. 



1995: IStone. Hawlev. Gammie. & Balbusl Il99f4 


Fleming & Stond 12003 


ISano. Inutsuka. Turner. & Stond 


2004: iGardiner & Stone 


2005). It is remarkable that this 



linear correlation between the stresses, in fully developed 
turbulent states resulting from very different sets of initial 
conditions, spans over six orders of magnitude. The ratio 
—Mrt/y/ ' Rr<t> = 4 is shown in the same figure with a dashed 
line, whereas the solid line shows the ratio obtained from 
equation 166^ for q = 3/2, i.e., —M r( p/R r( p = 5/3, charac- 
terizing the linear phase of the instability. Hence, to within 
factors of order unity, the ratio between the stresses in the 
turbulent state seems to be independent of the initial set 
of conditions over a wide range of parameter space and to 
be similar to the value set during the linear phase of the 
instability. 

The dependence of the ratio of the stresses on the shear 
parameter, q, has not been studied extensively with nu- 
meri cal simulations so far. The only co mprehensive study 
is by Hawlev, Balbus, & Winters (1999) and their result is 
shown in Figure 2] 3 . Superimposed on the figure is the an- 
alytic prediction, equation 16611 . for the ratio of the stresses 
as a function of the shear q at late times during the ex- 
ponential growth of the MRI. In this case, the qualitative 
trend followed by the ratio of the stresses at saturation as a 
function of the shear parameter q seems also to be similar 




[ , , . 1 

0.0 0.5 1.0 1.5 2.0 

q 

Figure 4. The dependence of the ratio of the mean Maxwell to 
the mean Reynolds stresses, — M r ^//? r 0, on the shear parameter 
q. The data points correspond to the results of shearing box sim- 
ulations by Hawlev. Balbus. fc Win ters! 1 19991) in the saturated 
state. The solid line shows the analytic result (eq. 1661 ) for the 
ratio of the mean stresses during the late time of the exponential 
growth of the MRI. 



to that obtained at late times during the linear phase of the 
instability. 

The simultaneous analysis of Figures [3] and 2] demon- 
strates that the ratio —M r $IR r $ during the turbulent sat- 
urated state in local simulations of accretion discs is de- 
termined almost entirely by the local shear and depends 
very weakly on the other properties of the flow or the ini- 
tial conditions. These figures also show that the ratios of 
the Maxwell to the Reynolds stresses calculated during the 
turbulent saturated state are qualitatively similar to the cor- 
responding ratios found during the late times of the linear 
phase of the instability, even though the latter are slightly 
lower (typically by a factor of 2). This is remarkable because, 
when deriving equation 1661 we have assumed that the MHD 
fluid is incompressible and considered only fluctuations that 
depend on the vertical, z, coordinate. Moreover, in the spirit 
of the linear analysis, we have not incorporated energy cas- 
cades between different scales, neither did we consider dissi- 
pation or reconnection processes that lead to saturation. Of 
course, the numerical simulations addressing the non-linear 
regime of the instability do not suffer from any of the ap- 
proximations invoked to solve for the temporal evolution of 
the stresses during the phase of exponential growth. Nev- 
ertheless, the ratio of the Maxwell to the Reynolds stresses 
that characterize the turbulent saturated state are similar 
(to within factors of order unity) to the ratios characterizing 
the late times of the linear phase of the instability. 



3 Notethat JlIawlev. Balbus. &: Winters! Jl999ft quote the average 
stresses and the width of their distribution throughout the simu- 
lations, but not the uncertainty in the mean values. In Fig.0] wc 
have assigned a nominal 30% uncertainty to their quoted mean 
values. This is comparable to the usual quoted uncertainty for 
the stresses and is also comparable to the spread in Fig. 151 
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APPENDIX A: MEAN VALUES OF 
CORRELATION FUNCTIONS 

For convenience, we provide here a brief demonstration 
of equation l|49|l . The relationship between the mean value of 
the product of two real functions, fg, and their correspond- 
ing Fourier transforms can be obtained by substituting the 
expressions for f(z,t) and g(z,t) in terms of their Fourier 
series, i.e., 



f(z,t) = /(*»,*) e* 



into the expression for the mean value, 



(Al) 



fg(t) 



2H 



f(z,t)g(t,z)dz. 



The result is, 



°° 1 



n.m.— — oo 



2H 



(A2) 



dz, (A3) 



which can be rewritten using the orthogonality of the Fourier 
polynomials in the interval [-H, H] as 



/<?(*) = E f(kn,t)g*(k n ,t). 



(A4) 



n— — oc 



This is the discrete version of Plancherel's theorem which 
states that the Fourier transform is an isometry, i.e., it 
prese rves the inner product (see., e.g., IShilov fc Silvermanl 

Denoting the Fourier transforms f(k n ,t) and g*(k n ,t) 
by f n and g* in order to simplify the notation, we can write 
the following series of identities 

oo 

fg(t) = E ( A5 ) 

71= — OO 

oo — 1 

= fo go + E f«-9n + Yl f n 9™ ( A6 ) 

n— 1 n — — oo 

oc 

= fo go + E + /-»$-»] (A7) 

n=l 

oo 

= fo go + [f n 9n + fnSn] (A8) 

n-l 

oo 

= MUgn], (A9) 

n = l 

where we have used the fact that the functions f(z,t) and 
g(z,t) are real and, hence, their Fourier transforms satisfy 
/_„ = and g- n = g^. Note that the factors fo and go are 
just the mean values of the functions f(z, t) and g(z, t) and, 
therefore, do not contribute to the final expression in equa- 
tion 14911 . It is then clear that, no matter whether the ini- 
tial Fourier transforms corresponding to the functions f(z, t) 
and g(z,t) are real or imaginary, the mean value, fg (t), will 
be well defined. 



APPENDIX B: ENERGETICS OF MRI-DRIVEN 
FLUCTUATIONS 

The relationships given in equation 1361 lead to identi- 
ties and inequalities involving the different mean stress com- 
ponents and the mean kinetic and magnetic energies asso- 
ciated with the fluctuations. In particular, the total mean 
stress is bounded by the total mean energy of the fluctua- 
tions. Moreover, the ratio of the mean magnetic to the mean 
kinetic energies is equal to the (absolute value of the) ratio 
between the mean Maxwell and the mean Reynolds stresses 
given by equation 166H . 

As we defined the mean stresses in terms of their per-k 
contributions in [21 we can a l so define, to leading order in 
time, the mean energies associated with the fluctuations in 
the velocity and magnetic field by 



E K (t) = 2j2£K(k n )\ ai \ 2 e 2 "<\ 



Em it) = 2^£ M (t) |ai| 2 e 



2-t 



(Bl) 



(B2) 



where the corresponding per-k contributions are given by 
E K (k„) = i[ftrr(ftn)+^(fcn)], (B3) 
£m{K) = ~ [Mrr(M +Mte(kn)]. (B4) 

Figure|2]shows the dependences of the functions £x(fcn) and 
£jw(fen) for a Keplerian profile, q = 3/2, and illustrates the 
fact that f Af(fcn) > £jf(fcn) for < k n < ksH and < q < 2. 

Using equation 13611 and the expression for the disper- 
sion relation (I2.'!li . it is easy to show that the following in- 
equalities hold for each wavenumber k n 

Hr&(k„) —Mr4>{k n ) 2~f(k n 



£i<(k n ) 



£M{k n ) 



^ 1. 



(B5) 



as long as q > 0. It immediately follows that the same in- 
equalities are also satisfied by the corresponding means, i.e., 



Rr^{t) ^ E K (t), 

-M r4> {t) <: E M (t). 



(B6) 
(B7) 



This result, in turn, implies that the total mean energy as- 
sociated with the fluctuations, E(t) = Ei((t) + Ej^it), sets 
an upper bound on the total mean stress, i.e., 



T r0 (t) < E(t) . 



(B8) 



At late times during the exponential growth of the instabil- 
ity, the growth of the fluctuations is dominated by the mode 
with k n = fc max and the mean stress T r( j,(t) will tend to the 
total mean energy E(t), i.e., 



livaTrJt) = \imE(t) . 



(B9) 



Furthermore, according to the first equality in equation 
1B5L we can conclude that the ratio of the mean magnetic 
to the mean kinetic energies has the same functional depen- 
dence on the shear parameter, q, as does the (negative of the) 
ratio between the mean Maxwell and the mean Reynolds 
stresses given by equation 1661 . i.e., 



E M it) 

i>>i E K {t) 



lim 



(BIO) 
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Therefore, the mean energy associated with magnetic fluctu- 
ations is always larger than the mean energy corresponding 
to kinetic fluctuations as long as the flow is Rayleigh-stable. 



REFERENCES 

Balbus S.A., 2003, ARA&A, 41, 555 
Balbus S.A., Hawley J.F., 1991, ApJ, 376, 214 
Balbus S.A., Hawley J.F., 1992, ApJ, 392, 662 
Balbus S.A., Hawley J.F., 1998, Rev. Mod. Phys., 70, 1 
Balbus S.A., Hawley J.F., 2002, ApJ, 573, 749 
Balbus S.A., Papaloizou J.C.B., 1999, ApJ, 521, 650 
Brandenburg A., Campbell C, Accretion Disks - New As- 
pects, Proceedings of the EARA Workshop Held in Garch- 
ing, Germany, 21-23 October 1996, XIII. Springer- Verlag, 
Berlin 

Brandenburg, A., Nordlund, A., Stein, R.F.; Torkelsson, 

U., 1995, ApJ, 446, 741 
Brandenburg A., 2001, ApJ, 550, 824 
Fleming T., Stone J.M., 2003, ApJ, 585, 908 
Gardiner T., Stone, J.M., 2005, in Magnetic Fields in the 

Universe, AIP Conf. Proc, 784, 475 
Goodman J., Xu G., 1994, ApJ, 432, 213 
Hawley J.F., Gammie C.F., Balbus S.A., 1995, ApJ, 440, 

742 

Hawley J.F., Gammie C.F., Balbus S.A., 1996, ApJ, 464, 
690 

Hawley J.F., Balbus S.A., Winters W.F., 1999, ApJ, 518, 
394 

Hawley J.F., 2000, ApJ, 528, 462 
Hawley J.F., 2001, ApJ, 554, 534 

Hoffman K.M., Kunze R., 1971, Linear Algebra, 2nd ed. 
Prentice Hall, N.J. 

Kato S., Fukue J., Mineshige S., 1998, Black-Hole Accre- 
tion Disks. Kyoto University Press., Kyoto 

Narayan R., Quataert E., Igumenshchev I.V., Abramowicz 
M.A., 2002, ApJ, 577, 295 

Sano T., Inutsuka S., Turner N.J., Stone J.M., 2004, ApJ, 
605, 321 

Shakura N.I., Sunyaev R.A., 1973, A&A, 24, 337 

Shilov G.E., Silverman R.A., 1973, Mathematical Analysis. 

MIT Press, Cambridge 
Stone J.M., Pringle J.E., 2001, MNRAS, 322, 461 
Stone J.M., Hawley J.F., Gammie C.F., Balbus S.A., 1996, 

ApJ, 463, 656 



© 2006 RAS, MNRAS 000, EHU 



